Optical and transport gaps in gated bilayer graphene 
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We discuss the effect of disorder on the band gap measured in bilayer graphene in optical and 
transport experiments. By calculating the optical conductivity and density of states using a micro- 
scopic model in the presence of disorder, we demonstrate that the gap associated with transport 
experiments is smaller than that associated with optical experiments. Intrinsic bilayer graphene has 
an optical conductivity in which the energy of the peaks associated with the interband transition are 
very robust against disorder and thus provide an estimate of the band gap. In contrast, extraction 
of the band gap from the optical conductivity of extrinsic bilayer graphene is almost impossible for 
significant levels of disorder due to the ambiguity of the transition peaks. The density of states 
contains an upper bound on the gap measured in transport experiments, and disorder has the effect 
of reducing this gap which explains why these experiments have so far been unable to replicate the 
large band gaps seen in optical measurements. 



Study of transport in graphene and related materials 
has recently become the focus of intensive research ef- 
forts since the experimental realization of gated samples 
in 2005 1 3 . Bilayer graphene (BLG) consists of two par- 
allel sheets of graphene coupled in the Bernal (AB) stack- 
ing arrangement, and the low-energy electronic quasipar- 
ticles behave like massive chiral fermions. This system 
shows unique physics, with many different physical ef- 
fects competing to dominate the observed properties and 
this has made BLG a very attractive material to study in 
the pursuit of understanding of fundamental physics of 
chiral materials. One of the most appealing features of 
this system is the dynamically tunable band gap which 
may be opened and controlled by a perpendicular electric 
field 4 '—. This is the most prominent reason for the explo- 
sion of interest in BLG because the potential applications 
of a tunable narrow gap semiconductor are legion. How- 
ever, the topic of charge transport through the gapped 
bilayer system is still a controversial issue. 

Various experiments have examined transport^ - — and 
the optical properties 1 ^"— of biased bilayer graphene. In 
many cases, the band gap extracted from a transport 
measurement was many times smaller than that found 
optically. The most obvious explanation for this dis- 
crepancy is sample disorder such as charged impurities 
and short-range scatterers 1 ^ 9 . Some theoretical work 
on the density of states (DOS) and optical conductivity 
of biased BLG in the presence of disorder using coherent 
potential approximation 20 , in the instanton approach to 
the in-gap fluctuation states^, for lattice defects^, for 
midgap states and Coulomb scatterers 2 ^, and for finite- 
ranged scatterers 24 has been published, but to this date 
no theoretical studies exist which consider the issue of 
transport and optical gaps in BLG in the presence of 
disorder. 

In this Rapid Communication, we present a compre- 
hensive theory of the optical conductivity in gapped BLG 
and the associated DOS. From these quantities we can de- 
scribe the gaps measured by optical and transport mea- 
surements and make direct comparisons between them 
for the same system. The optical gap is best found from 



intrinsic graphene (where there are no excess charge car- 
riers) because the relevant features in the optical con- 
ductivity are much more robust against disorder in that 
case. In contrast, the measurement of the band gap in 
extrinsic graphene (i.e. BLG with finite carrier density) 
relies on spectral features which are strongly affected by 
disorder. One of the main achievements of this work is 
the development of a quantitative microscopic theory for 
the optical conductivity of BLG in the presence of dis- 
order, and the direct comparison with the DOS at the 
same level of approximation. We also demonstrate for 
the first time why the gaps extracted via different ex- 
perimental techniques vary so widely. Finally, we shall 
propose that manufacturing very clean samples will lead 
to much closer agreement in the optical and transport 
gaps. 

To accomplish this, we compute the DOS and optical 
conductivity including disorder within the self-consistent 
Born approximation (SCBA). The results will be sim- 
ilar for various scattering mechanisms, but we choose 
Coulomb scatterers with Thomas- Fermi screening 25 ^ to 
illustrate them. Therefore, the impurity density is a pa- 
rameter which describes the effective disorder concentra- 
tion, and not strictly the density of charged impurities. 

In the SCBA, the self-energy is computed as follows: 
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Here A is the band index, ni mp is the impurity density, r] 
is an infinitesimally small positive number, F\A'(k, k') 
is square of the wave function overlap between |A,k) 
and |A',k') states, which is a nontrivial function of 
the wave vectors and band indices. The Fourier trans- 
form of the impurity potential is denoted by Vi mp (q) at 
wave vector q. This self-energy can be used to com- 
pute the retarded electron Green's function G\(\z, E) = 
(E - £ A k - S A (k, E) + if])' 1 in the presence of disorder, 
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from which the DOS can be extracted in the usual way: 
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This includes the spin and valley degeneracy factors g s = 
g v = 2 since the charged impurity scattering does not 
mix valleys or spins. 

For the screened Coulomb interaction with the effective 
impurity distance di mp from the graphene layer, VJ mp (q) 
is given by 
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where eo is the background dielectric constant and q s is 
the screening wave vector. For a finite carrier density 
(the extrinsic case), the screening wave vector can be 
approximated by q s ~ <2tf where #tf = ^--DoC^f) is 
the Thomas- Fermi wave vector and Dq(E-f) is DOS at the 
Fermi energy in the absence of disorder. However, 
if the Fermi energy lies in the gapped region (i.e. for 
intrinsic graphene), there is no DOS at the Fermi energy 
and the screening wave vector should be defined by the 
interband contribution to the polarization. These two 
methods of defining the screening wave vectors provide 
similar screening lengths 27 . 

The optical conductivity is computed within the Kubo 
formalism 2 ^ and at zero temperature can be summarized 
by the following expression: 
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and (j) is the angle of the wave vector k. Throughout 
this Rapid Communication, we use the four band tight- 
binding Hamiltonian 3 -^ of BLG which includes the in- 
tralayer nearest-neighbor hopping elements parametrized 
by 7o = 3 eV, the hopping via interlayer dimer bonds 
7i = 0.3 eV, and the onsite potential energy supplied 
by the gates given by ±u/2 depending on which layer 
contains the lattice site 3 . 

We can now describe the transport and optical gaps 
within the framework of our theory, defining the gaps as 
follows. We characterize the gap in the absence of disor- 
der by the associated potential energy difference between 
the layers in the single particle tight-binding formalism, 
and label this quantity by u. The "sombrero" dispersion 
of gapped BLG 4 - means that the minimum band gap oc- 
curs slightly away from the K point such that the bare 
band gap is E^ p = ^7i/V7i + u 2 . The DOS gap is a 




FIG. 1. (Color online) Density of states (DOS) of intrinsic bi- 
ased BLG with u = 0.2 eV on Si02 substrate. The case with 
no charged impurities is shown by the black line. The other 
lines show the DOS for increasing charged impurity concen- 
tration, with n imp =l, 5, 10, 20xl0 12 cm~ 2 . 



theoretical quantity defined as the energy difference be- 
tween the valence and conduction bands as defined by 
the onset of the DOS in the SCBA calculation (see the 
inset to Fig. [1]), while the transport gap is extracted from 
experiments. We note that, due to the possible existence 
of mid-gap impurity states as well as puddle formation 29 , 
the transport gap may be lower than the DOS gap. To 
the best of our knowledge, no evidence of strong local- 
ization has so far been seen in transport experiments on 
BLG, and this allows us to assume that the DOS gap is 
the effective upper bound for the transport gap. 

Following the experimental definition 14 , the optical 
gap is defined as the energy at the top of the first peak in 
the optical conductivity of intrinsic graphene. This peak 
results from the direct optical excitations from the va- 
lence to the conduction band and is the most direct way 
of extracting the gap from an optical experiment. This 
gap is illustrated in the inset to Fig. [3fb). 

The effect of disorder on the DOS gap is easy to under- 
stand in Fig. [TJ With increasing impurity concentration, 
the gap size reduces and eventually closes at a finite value 
of the impurity concentration. This figure shows the DOS 
for a fairly large potential asymmetry u = 0.2 eV (which 

gives £?gap = 0.166 eV), and a smaller gap generated by 
a smaller asymmetry would be closed by a smaller con- 
centration of charged impurities. 

The effect of disorder on the optical gap is also easy 
to determine if it is extracted from the optical conduc- 
tivity of intrinsic BLG. Figure [2fa) shows the interband 
contribution to the nondisordered optical conductivity of 
intrinsic BLG 30 along with a sketch showing the origin of 
the two peaks. The transition labeled "0" is the one with 
the lowest energy and corresponds to the excitation of 
carriers across the band gap. Figure [2fb) shows the same 
quantity in the presence of disorder. In addition to in- 
cluding the impurity scattering effects from ni mp ^ 0, we 
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FIG. 2. (Color online) (a) Interband optical conductivity 
of biased intrinsic BLG with u — 0.2 eV in the absence of 
disorder. The inset shows the band structure with the in- 
terband transitions indicated by arrows, (b) Optical con- 
ductivity of disordered biased BLG for impurity densities 
n imp 0, 1, 5, 10, 20 x 10 12 cm" 2 . 
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FIG. 3. (Color online) DOS gap and optical gap in intrin- 
sic BLG as a function of impurity density. Parameters in 
(a) correspond to an Si02 substrate and in (b) to a h-BN 
substrate^. To aid comparison, the charged impurity con- 
centration of an Si02 substrate is of the order of 10 12 cm -2 
whereas for a h-BN substrate it is of the order of 10 11 cm -2 . 



have included a finite phenomenological level-broadening 
of j] = 1 meV in the numerics which slightly rounds 
off the peaks in the ni mp = case. The introduction 
of charged impurities further broadens the peaks corre- 
sponding to the interband transitions but leaves their 
position almost unchanged. This explains why the opti- 
cal gap extracted from experiment [which is defined by 
the peak value of ct(uj)] is very close to that predicted by 
tight-binding theory without disorder. 

The DOS and optical gaps extracted from these cal- 
culations are shown in Fig. [3] as a function of the im- 
purity density 31 . The gaps are identical in the nondis- 
ordered system, but as the density of charged impurities 
increases, the DOS gap is suppressed much more quickly 
than the optical gap. In the majority of the original ex- 
periments on transport in BLG, the samples are made 
with exfoliated graphene on Si02 substrates and would 
therefore have had a large amount of charged impurity 
disorder. We believe that this observation accounts for 
the difference in the measured optical and transport gaps 



FIG. 4. (Color online) (a) Interband optical conductivity of 
extrinsic biased BLG with carrier density n = 5 x 10 12 cm -2 
and u = 0.2 eV in the absence of disorder. The inset shows 
the band structure with the interband transitions indicated by 
arrows and the dotted line shows the Fermi energy, (b) Opti- 
cal conductivity of disordered biased BLG for several values 
of the impurity density. 



in this system. The actual transport gap could be con- 
siderably smaller than our theoretical DOS gap 29 , as we 
have already noted. 

Figure [H shows the optical conductivity calculated for 
extrinsic BLG on an Si02 substrate 3 ^. In this case, 
the extra transitions [which are sketched in the inset to 
Fig. IK a)] make the optical conductivity a highly nontriv- 
ial function of energy and no simple way of extracting 
an effective optical gap seems obvious. The interband 
optical conductivity (calculated using the formula from 
Ref. l30l ) in the nondisordered case is shown in Fig. IHa). 
The inclusion of disorder [FigSJ^b), using Eq. (jH] further 
masks the structure and makes it very difficult to reliably 
identify the peaks in the observed optical spectra. The 
peak at low energy is due to the Drude contribution from 
intraband processes near the Fermi surface. For exam- 
ple, several experimental studies extract the size of the 
band gap by comparing the energy of two peaks relating 
to different interband transitions. In the language of Fig. 
[H(a), these transitions are labeled 1 (purple online) and 
5 (red online). Making the identification with the plot 
in Fig. HJb), we see that transition 1 corresponds to the 
second (smaller) peak, while transition 5 corresponds to 
the slight shoulder on the low energy side of the main 
peak. As disorder increases, the peak due to transition 
1 quickly disappears, and the shoulder due to transition 
5 becomes indistinct. It is therefore impossible to use 
this method to extract the optical gap in the presence of 
a realistic amount of charged impurity disorder in BLG, 
and the peaks seen in the optical conductivity are not 
obviously related to any individual transition. This is a 
qualitative new finding of our work. 

Finally, we discuss a practical way to test our theory 
and to bring the DOS and optical gaps much closer in en- 
ergy. To do this, optical spectroscopy and charge trans- 
port experiments should be performed on the same BLG 
device using an ultra-clean substrate, such as hexagonal 
boron nitride (h-BN) which has an atomically smooth 
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FIG. 5. (Color online) Optical conductivity of biased BLG 
with u — 0.2 eV on different substrates, (a) ni mp = 5 x 
10 12 cm -2 on Si02 substrate and (b) ni mp = 10 11 cm -2 on an 
h-BN substrate. In each case, the black line shows the intrin- 
sic case, and the other lines show different carrier densities 
for extrinsic graphene, measured in 10 12 cm -2 . 

surface and much lower impurity concentration than cur- 
rent Si02 device a 32 i 33 . Recently, gated graphene layers 
on h-BN substrates have been successfully fabricated 32 , 
and it was demonstrated that graphene on h-BN sub- 
strates have substantially higher mobility than graphene 
on Si02 substrates by an order of magnitude or more. 

Figures [SJa) and [5jb) show the optical conductivity 
for a typical impurity density of Si02 substrate and h- 
BN substrate for a range of carrier densities. Apart from 
different disorder strengths, the two substrates are distin- 
guished by different dielectric environments 3 ^. Because 
of the relatively low impurity density, the optical conduc- 
tivity for the h-BN substrate shows sharper peaks than 



those of Si02 substrate allowing for more accurate iden- 
tification of the origin of the various features and hence 
for a more accurate determination of the band gap. Fig- 
ure [31(b) shows energy gaps extracted from the DOS and 
optical conductivity of biased intrinsic BLG with u = 0.2 
eV on an h-BN substrate as a function of impurity den- 
sity range. In contrast to the same plots for Si02 shown 
in Fig. [3](a), the deviation of the DOS gap from the opti- 
cal gap for the h-BN substrate is rather small. Note that 
experimentally realistic values of the impurity density on 
the h-BN substrate 33 are of the order of 10 11 cm -2 so 
that the DOS gap is approximately 95% of the optical 
gap. Therefore, BLG mounted on an h-BN substrate can 
be used as a tunable gap semiconductor and will yield a 
significantly higher on/off ratio than the corresponding 
Si02 substrate. 

In summary, using the SCBA, we have shown that op- 
tical experiments and transport measurements do not al- 
ways probe the same gap in biased BLG. The disorder 
affects the DOS so that the gap is reduced, while the po- 
sition of the peak associated with interband transitions in 
the optical conductivity of intrinsic BLG is rather robust 
against disorder. In contrast, the features in the opti- 
cal conductivity of extrinsic BLG which are compared 
to extract the band gap become indistinct with signifi- 
cant disorder. Our work clearly establishes the danger of 
extracting a na'ive optical band gap in bilayer graphene 
from the optical absorption data without a thorough the- 
oretical analysis, calling into question several recent ex- 
perimental claims^ 3 -"— . We also suggest that samples 
with low disorder (such as BLG on an h-BN substrate, 
or suspended BLG) should show comparable DOS and 
optical gaps. 
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